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Abstract. The VERITAS collaboration reports the detection of very-high-energy (VHE) gamma-ray emission from the high- 
frequency-peaked BL Lac object lES 1218+304 located at a redshift of z = 0.182. A gamma-ray signal was detected with 
high statistical significance for the observations taken during several months in the 2006-2007 observing season. The photon 
spectrum between ~ 160 GeV and ~ 1.8 TeV is well described by a power law with an index of F = 3.08 ± 0.34s/ „/ ±0.2_5yj. 
The integral flux above 200 GeV corresponds to --^ 6% of that of the Crab Nebula. The light curve does not show any evidence 
for VHE flux variability. Using lower limits on the density of the extragalactic background light (EBL) in the near-IR to mid- 
IR we are able to limit the range of intrinsic energy spectra for lES 1218+304. We show that the intrinsic photon spectrum is 
harder than a power law with an index of F = 2.32 ± 0.37. When including constraints from the spectra of lES 1 101-232 and 
lES 0229+200, the spectrum of lES 1218+304 is likely to be harder than F = 1.86 ±0.37. 
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INTRODUCTION 

In the TeV energy range, 19 BL Lac objects and 1 Ra- 
dio Galaxy (M 87) have been established as emitters of 
TeV gamma rays. With the exception of the recently dis- 
covered low-frequency-peaked BL Lac (LBL) object BL 
Lacertae Ol and the intermediate-frequency-peaked BL 
Lac (IBL) object W Comae [iZ], all TeV blazars are high- 
frequency-peaked BL Lac (HBL) objects. HBLs are a 
subclass of blazars characterized by a synchrotron peak 
at X-ray energies, unlike quasars that have higher lumi- 
nosity and a synchrotron peak at optical/IR energies. lES 
1218+304 is an X-ray -bright (FikeV > 2^Jy) HBL ob- 
ject located at a redshift z — 0.182 fs"]. Based on SED 
modeling and BeppoSAX X-ray spectra, several of the 
HBLs were predicted to be TeV sources and several of 
them have indeed been detected at TeV energies f?]. lES 
1218+304 was predicted to be a good TeV candidate 
based on the position of its synchrotron peak at high en- 
ergy and sufficient radio-to-optical flux. It was recently 
detected (at very high energies) by the MAGIC telescope 
11. 



VERITAS OBSERVATIONS 

The VERITAS observatory consists of an array of four 
12-meter diameter imaging atmospheric Cherenkov tele- 
scopes (lACTs) located at the Fred Lawrence Whipple 
Observatory in southern Arizona |@]. The array is sensi- 
tive to gamma rays in the energy range from 100 GeV 
to 30 TeV. The first two telescopes were operated in 
the stereoscopic observation mode from March 2006 and 
the third and fourth telescopes came online in Decem- 



ber 2006 and April 2007, respectively VERITAS ob- 
served lES 1218+304 from December 2006 to March 
2007 using three telescopes. The data were taken in wob- 
ble mode where the source is offset from the center of the 
field of view by 0.5° and the background is measured di- 
rectly from different regions in the same field of view but 
away from the source region. After removing data taken 
under poor sky conditions or affected by various detector 
problems, we were left with a total observation time of 
17.4 hours covering a range in zenith angle from 2° to 
35°, with an average zenith angle of 14°. 

ANALYSIS AND RESULTS 

The standard analysis tools were used to process the data. 
After calculating the standard Hillas parameters, images 
with an integrated charge less than 400 digital count^ 
or with a distance from the center of the camera larger 
than 1 .2° were rejected. The shower direction in the field 
of view and the impact parameter of the shower core 
were calculated using stereoscopic techniques ||7|, ISj]. 
The background was estimated using the reflected-reeion 
and ring background methods as described in j^. A 
set of scaled cuts on the width and length parameters 
were used to identify gamma-ray events in the data. An 
extensive set of Monte Carlo simulations was used to 
generate lookup tables and calculate the effective area of 
the detector as a function of zenith angle and gamma-ray 
energy. 



An integrated charge of 400 digital counts calculated using a 10 ns 
integration window corresponds to ~75 photoelectrons. 
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FIGURE 1. Significance map of tlie region around lES 
1218+304. The white cross indicates the position of the radio 
source. The white circle in the upper left corner shows the angu- 
lar resolution of VERITAS (PSF). The ring background model 
was used to estimate the background and the significances were 
calculated using equation 17 in Li and Ma 1. 10.1 . 

Figure [1] shows a two-dimensional sky map of signif- 
icances for a region centered on the radio coordinates of 
lES 1218H-304 (12''21'"21.9',+30°10'37" J2000). The 
excess is compatible with that expected from a point 
source. Fitting a two-dimensional Gaussian distribution 
to the uncorrelated excess map yields a peak position in 
good agreement with the radio coordinates. 

Figure |2] shows the light curve of the integral flux 
above 200 GeV for the months of January, Febru- 
ary, and March 2007. The average integral flux was 
calculated for each day assuming a spectral shape of 
dN/dE oc with T = 3.08. A statistical test (x^/dof 
= 6.7/11) indicates that no statistically significant vari- 
ability was detected in the data. Figure [3] shows the 
time-average differential energy spectrum for gamma- 
ray energies between 160 GeV and 1.8 TeV. The 
shape is consistent with a power law (x^/doi = 2.1/5) 
dN/dE = C X (E/0.5TeV)"'^ with a photon index T = 
3.08 ± 0.34s,ar ± 0.2jyf and a flux normalization con- 
stant C = (7.5 ± 1.1,, „, ± 1.5,,.,) X 10-'2cm-2s-iTeV-^ 
The integral flux is <I>(E > 200 GeV) = (12.2 ± 2.6) x 
lO^^^cm^^s^' which corresponds to ^ 6% of the flux of 
the Crab Nebula above the same threshold. 



DISCUSSION AND CONCLUSIONS 

Due to the relatively large redshift of lES 1218H-304 (z = 
0.182) we expect significant attenuation from the inter- 
action of high-energy gamma rays with low-energy pho- 
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FIGURE 2. Light curve of the integral photon flux above 200 
GeV for the source lES 1218+304. Each point corresponds to 
the average daily flux and assumes a spectral shape of dN/dE 
oc E^'" with r = 3.08. The error bars represent the statistical 
uncertainty. 



tons from the EBL. TeV gamma-rays from extragalac- 
tic sources are expected to pair produce with optical/IR 
low-energy photons as they travel through intergalactic 
space, leading to a cutoff in the measured gamma-ray 
spectrum. Intrinsic spectra of blazars are not known a 
priori. Nevertheless, it is possible to use the EBL lower 
limits from galaxy counts to restrict the range of EBL 
scenarios / models and discern a coiTesponding range of 
intrinsic blazar spectra. 

Here we present possible intrinsic energy spectra of 
lES 1218+304 that are compatible with EBL lower hm- 
its from galaxy counts. Furthermore, we broaden our 
study and include the energy spectra of lES 1101-232 
and lES 0229+200 as measured by H.E.S.S. (see UHl^ 
respectively) to reach a conclusion about the intrinsic 
energy spectra of this small sample of TeV blazars. 
The combination of energy spectra in the sub-TeV to 
TeV waveband at large redshfift (lES 1102-232, lES 
1218+304) with energy spectra in the multi-TeV regime 
at moderate redshift (lES 0229+200) provides additional 
sensitivity to EBL spectra and their relative intensities in 
the near-IR and mid-IR. We have considerd a large vari- 
ety of EBL scenarios with different spectral indices and 
shapes. A convenient parameterization of the EBL sce- 
narios, provided by il4ll . is used in this study for provid- 
ing a limit to the hardness of the blazar spectrum of lES 
1218+304 as it provides a wide range of EBL spectra 
with different near-IR to mid-IR ratios. 

Table [T] shows the intrinsic energy spectra of all 3 
blazars for a range of EBL scenarios. The ones that are 
still compatible with lower limits from galaxy counts 
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FIGURE 3. Differential energy spectrum of VHE photons 
above 160 GeV for lES 1218+304. The markers indicate 
measured data points and the continuous line is a power-law 
fit. Downward pointing arrows correspond to upper flux lim- 
its (99% probability, Helene |11]) for bins with significances 
below two standard deviations. The dot-dash line shows the 
differential energy spectrum of lES 1218+304 measured by 
MAGIC 1 5] and the dash line shows the differential energy 
spectrum of the Crab Nebula measured by VERITAS (Septem- 
ber to November 2006) for comparison. 



represent the softest possible gamma-ray spectra. When 
considering lES 1218+304 by itself, the softest intrinsic 
spectrum is described by a power law with F = 2.32 ± 
0.37 and is derived from a scaled version of the LHL 
scenario (LHL0.76). However, when applying this anal- 
ysis to previously reported blazar spectra lES 1 101-232 
and lES 0229+200, the LHL0.76 scenario would require 
an extremely hard intrinsic spectrum for IBS 0229+200. 
Therefore, in order to provide a limit to the hardness of 
the blazar spectra based on these three sources we search 
for the softest possible intrinsic spectrum. As can be seen 
in Table [H a search for the softest possible blazar spec- 
trum among the sample of three blazars yield an EBL 
scenario (AHA x 0.45) that still requires the spectrum of 
lES 1218+304 to be as hard as r= 1.86 ±0.37. All other 
EBL scenarios yield harder spectra for one of these three 
blazars. 

In conclusion, these results indicate that blazar spectra 
as evidenced by the two most distant blazars, lES 1101- 
232 and lES 1218+304, are hard. These values are still 
within the acceptable range predicted for shock acceler- 
ation in blazars (see, e.g. II15I1 ). however they are also 
close to a limit (F — 1.5) that was previously suggested 
by lfT2ll . The results presented here from VERITAS ob- 
servations confirm with high statistical significance the 
MAGIC discovery of the HBL object lES 1218+304 
as a source of VHE gamma rays. The detection of 



lES 1218+304 together with the other blazars studied by 
VERITAS lfl6ll give valuable insights into the physics of 
these sources and the EBL. 
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TABLE 1. The absorption-corrected spectral indices (dN/dE 
E-rso„roc) of lES 1101-232, lES 1218-^304, and lES 0229-1-200 from 
a small sample of all EBL scenarios considered. EBL scenarios that 
fall below the lower limits from galaxy counts are marked with an 
asterisk. Scaled Aharonian scenarios (AHA.i:.j) are taken from 1 12]. 
HHH (High near-IR, mid-IR, far-IR)m MHL (Medium near-IR, High 
mid-IR, Low far-IR), etc, are taken from |14]. 



Scenario 




^source 




lES 1101-232 


lES 1218+304 


lES 0229-1-200 


AHA 1.0 


0.44 ±0.20 


0.35 ±0.40 


2.35±0.14 


AHA0.55 


1.54±0.20 


1.59±0.37 


2.41 ±0.14 


AHA0.45 


1.78±0.20 


1.86±0.37 


2.43±0.13 


HHH 


-0.67 ±0.12 


-0.72 ±0.29 


0.90±0.17 


LLH* 


2.01 ±0.22 


2.07 ±0.35 


2. 12 ±0.20 


LHL 


2.04 ±0.20 


2.08 ±0.39 


0.94 ±0.32 


LHL0.70* 


2.32±0.21 


2.43 ±0.37 


1.43 ±0.29 


LHL0.76 


2.23 ±0.21 


2.32 ±0.37 


1.30±0.29 


LHL0.82 


2.18±0.21 


2.26 ±0.38 


1.20±0.30 


MHL0.70 


1.26±0.19 


1.34±0.36 


1.35±0.21 


MHL0.55 


1.61±0.19 


1.73 ±0.35 


1.59±0.20 


LLL* 


2.06±0.16 


2.20±0.34 


2.11 ±0.20 



